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Abstract 

The configuration-fixed deformation constrained relativistic mean field approach with time-odd 



(N 



component has been applied to investigate the ground-state properties of 33 Mg with effective 
interaction PK1. The ground state of 33 Mg has been found to be prolate deformed, = 0.23, with 
the odd neutron in 1/2 [330] orbital and the energy —251.85 MeV which is close to the data —252.06 



> 
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MeV. The magnetic moment — 0.9134 un is obtained with the effective electromagnetic current 



which well reproduces the data — 0.7456 self-consistently without introducing any parameter. 
The energy splittings of time reversal conjugate states, the neutron current, the energy contribution 
from the nuclear magnetic potential, and the effect of core polarization are discussed in detail. 
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I. INTRODUCTION 



Nuclear magnetic moment, as one of the most important physical observables, provides 
key information to understand nuclear structure, and has attracted the attention of nuclear 
prists since the ea rly da y8 M- Ap-t f r o,n tee mag „et,e m0m eet S of 8tab ,e nuclei Q, 
it is now even possible to measure the nuclear magnetic moments of many short-lived nuclei 
far from the stability line with high precision [s| with the development of the radioactive 
ion beam (RIB) technique. 

Recently the nuclear magnetic moments of 33 Mg has become a hot topic due to the 
following reasons: 1) it is a neutron- rich nucleus close to the so-called "island of inversion" f| ; 
2) different spins and configurations for the ground state of 33 Mg are assigned in a series of 



experiments 



Sms an 
lflB n 



101 ]. In order to remove the confusion, the spin and magnetic moment 
for the ground state in 33 Mg have been directly measured in Ref. ll| with I = 3/2 and 
\i = — 0.7456(5)yUN, which becomes a test for various theoretical approaches. In shell-model, 
the magnetic moment of the ground-state in 33 Mg, can be reproduced only in the model 
space with 2p — 2h confi gur ation ll| . With the assignment of configuration in Ref. lOj, the 
simple Additivity Rules 12J can only account for half of the experimental magnetic moment. 

For the last two decades, the relativistic mean field (RMF) theory 13J has achieved 
great success in describing many phenomena in the nuclei near and far from the line of 

/3-stability HQ" 



16| . In the widely used version of RMF models, only the time-even 



component of vector meson fields are taken into account. In odd-A or odd-odd nuclei, 
however, the baryon current due to the unpaired valence nucleons will lead to the time-odd 
components of vector fields, i.e., the time-odd fields. The time-odd fields together with the 
corresponding core polarization will modify the nuclear current, single-particle Dirac spinor, 



and magnetic moments, etc. The importance of the time-ode 



nc 



13, 



ds has been demonstrated 



19| , moment of inertia for 



in the successful descriptions of nuclear magnetic moments 

identical superdeformed bands [^J, and Ml transition rates in magnetic rotation 2ll. l22|. 
etc. 

In this paper, the ground-state properties of 33 Mg including the binding energy and 
magnetic moment will be investigated in axially deformed RMF approach with time-odd 



components of vector meson field 



In order to examine the magnetic moments for 



different valence nucleon configurations, both adiabatic and configuration-fixed deformation 



2 



constrained calculation 



24j will be performed. 



II. THEORETICAL FRAMEWORK 



The basic ansatz of RMF theory is a lagrangian density jl3, 14, If], |25| where nucleons are 
described as Dirac particles which interact via the exchange of various mesons and photon. 
The lagrangian density is written in the form 

-\g^ - in^O^ + im>^ + ^c 3 (oA, M ) 2 

■ R»» + ■ h - \f" v f^ (i) 

The meson fields contain the isoscalar a meson, isoscalar-vector u meson and the isovector- 
vector p meson. M and rrii(gi) (i = a,cu,p) are the masses (coupling constants) of the 
nucleon and the mesons respectively and 

W = (2a) 
RT = Vp-d"? (2b) 
= d»A u -d u A» (2c) 

are the field tensors of the vector mesons and the electromagnetic field. In this paper, we 
adopt the arrows to indicate vectors in isospin space and bold type for the space vectors. 

For even-even nuclei with time reversal symmetry, the time-odd components of vector 
mesons and photon fields do not contribute to the energy functional. In odd-A nuclei, the 
odd nucleon breaks the time-reversal invariance, and time-odd fields give rise to the nuclear 
magnetic potential, i.e., the time-odd component of the vector potential. Then the equation 
of motion for nucleon can be obtained as: 

{a • [-zV - V(r)} + V (r) + P[M + S(r)]}ifj t = e^, (3) 

with the scalar potential S(r) = g a a(r), the time-like component of vector potential Vo(r) = 
g w u) (r) + g p T 3 p (r) + e^^-Ao^r), and the time-odd component of vector potential V(r) = 
p w a>(r), where p{r) and A(r) are neglected since they turn out to be small compared with 



u?(r) field in light nuclei 17]. 



The Klein-Gordon equations for a, time-like components of vector mesons fields ujq, po 



16|. The time-odd component of 



and electromagnetic fields Aq are the same as in the Ref 
uj meson is determined by 

{-A + ml }uj = g UJ j B - c 3 u"u v u, (4) 
with the baryon current j B = ^n^-y^. The summation is confined to the particle states 

i 

with positive energies in the no-sea approximation. As the pair correlation is neglected here, 
the occupation numbers rii take the value one (zero) for the states below (above) the Fermi 
surface. Restricting to an axially symmetric representation here, only azimuthal baryon 
currents j B (z, r±) on circular lines around the symmetry . , only the 

azimuthal component of the uj vector field exists. 

The total energy of the system, including the time-odd fields, is, 

E = -Epart + E a + Eu + E p + E c + E c , m , (5) 

with 

Eu = -7; J d 3 r \g u ji(T)u (r) - ^^(r)j'^(r) 

-^c 3 K^) 2 (r)] (6) 



and energy contributions from other parts, not shown here, are the same as given in Ref. 



19]. 



To describe the nuc 
which is defined as 118 



ear magnetic moment, one needs the effective electromagnetic current 



>(x) = ^(x^l^^x) + ^[^(*K>(x)], (7) 

where a^ v = |[7 M ,7 y ], and K is the free anomalous gyromagnetic ratio of the nucleon: 
k v = 1.793 and k" = —1.913. The spatial component of the effective electromagnetic 
current operator is given by 

j(r) = ^+( r) al_^^(r) + x [V>+(r)0Ity(r)], (8) 

where the first term is the Dirac current, and the second term is the so-called anomalous 
current. 

From the effective electromagnetic current, the magnetic moment can be obtained by 

H = - J d 3 rrx (g.s.|j(r)|g.s.) . (9) 
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III. NUMERICAL DETAILS 



The Dirac equation for nucleons and the KleinCGordon equations for mesons are solved 
by expansion in terms of the isotropic harmonic oscillator basis functions in cylindrical 
coordinates with 14 oscillator shells for both the fermion and the boson fields 26|, |27j. The 
oscillator frequency is given by hwo = 41A -1 / 3 MeV. In odd A nuclei, as the time reversal 
invariance is broken by the odd nucleon, the Dirac equation for nucleons should be solved 
separately in two subspaces {ipj} and {ipj}, where ipj and tpj are time reversal conjugate 



states 



19|]. 



The energy surface is obtained throu gh t he deformation constrained calculation. The 



quadratic constraint is adopted as in Ref. 
(Q2) to a given value q, i.e., 



281 ] by constraining the mass quadrupole moment 



1 



(H') = (H) + -C((Q 2 )-q)\ 



(10) 



where (H) is the total energy, and C is the stiffness constant. The deformation parameter 
P2 is obtained from the calculated (Q 2 ) for the protons and neutrons through 

3 



(Q2) = <Q 2P > + (Q*n) 



(11) 



with Rq = 1.2 A 1 / 3 . In the following, both the adiabatic and the configuration-fixed defor- 



mation constrained calculation 



24 



301 ] will be p erformed. 



The effective interaction parameter set PK1 3j| is used throughout the calculation and 
the center-of-mass (cm.) correction is taken into account microscopically by 

1 



El 



2MA 



(P 2 ) 

\ x cm./ i 



(12) 



where P c . m . is the total momentum operator of a nucleus with mass number A. 



IV. RESULTS AND DISCUSSION 

The energy surfaces for 33 Mg as a function of the quadrupole deformation parameter fa 
calculated by adiabatic, shown as open circles, and configuration-fixed, shown as solid lines, 
deformation constrained RMF approach with time-odd component using PK1 parameter set 
are presented in Fig. Ufa). As discussed in the Ref. 3(]], the configuration-fixed deformation 
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FIG. 1: (Color online) (a) The energy surfaces for 33 Mg as a function of 02 by adiabatic (open 
circles) and configuration-fixed (solid lines) deformation constrained RMF approach with time-odd 
component using PK1 parameter set. The minima in the energy surfaces for fixed configurations are 
represented as stars and respectively, labeled as A, B, C, D, E, F, and G of which A is the ground 
state with total energy E = —251.85 MeV and 02 = 0.23, in comparison with the data —252.06 
MeV [32| (dotted line), (b) Magnetic moments for the corresponding configurations in panel (a) 



as a function of 02- The magnetic moment for the ground state is — 0.913 (in in comparison with 



the experimental value fJ, = — 0.7456(5) /xjv 



ll|] (dotted line). 



constrained calculation gives a continuous and smooth curve for the energy surfaces as a 
function of 02- The local minima in the energy surfaces for each configuration are represented 
by stars and labeled as ACG in ascending order of energy. The ground state A is found to be 
prolate deformed, 02 = 0.23, with the total energy —251.85 MeV, which is close to the data 
—252.06 MeV 32|. Using Eq. (Q, the effective electromagnetic current gives the nuclear 
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FIG. 2: (Color online) Neutron single-particle energies for 33 Mg as a function of obtained by 
configuration-fixed deformation constrained calculation for the configuration of ground state A. 
Positive (negative) parity states are marked by solid (dashed) lines. Each pair of time reversal 
conjugate states splits up into two levels with the third component of total angular momentum 
fi > and Cl < denoted by red and black lines respectively. The solid circle denotes that the 
corresponding orbitals are occupied in the ground state. 

magnetic moment for given configurations in Fig.QTb). It is found that the magnetic moment 
is sensitive to the configuration, but not so much to The magnetic moment for the ground 



state is — 0.913 hn which is in good agreement with the data \x = — 0.7456(5) \x n 111 ] . 



compared with the shell-model results — 0.675 /xjv and — 0.705 [an U| restricted to 2p — 2h 
configuration using two different interactions designed specifically for the island of inversion. 

In order to examine the evolution of the single particle level and compare with the re- 
sults in Ref. Ill ], neutron single-particle energies for 33 Mg as a function of (32 obtained 
by configuration-fixed deformation constrained calculation for the configuration of ground 
state A are presented in Fig. [2] The positive (negative) parity states are marked by solid 
(dashed) lines, and the occupied orbitals are represented by filled circles. The self-consistent 
calculation here gives the odd neutron in 1/2 [330] orbital with /3 2 = 0.23 for the ground 



state, while in Ref. 11|, in a Nilsson-model picture, the odd neutron in 3/2[321] orbital with 



prolate deformation 0.3 < (3% < 0.5 is proposed to reproduce the spin and parity I n — | . 

As the time reversal invariance is broken by the nuclear magnetic potential, each pair of 
time reversal conjugate states splits up into two levels with the third component of total 
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FIG. 3: Azimuthal components of neutron current in the xy plane at z = 0.40 fm (a), z = 1.20 
fm (b), z = 2.0 fm (c), z = 2.81 fm (d), z = 3.63 fm (e), and z = 4.48 fm (f) respectively for the 
ground state of 33 Mg. The direction and length of the arrows respectively represent the orientation 
and magnitude of current which is in different scale for different panels. 

angular momentum Q > and Cl < 0. The energy splittings of time reversal conjugate 
states range from 0.01 to 0.2 MeV, and the largest splitting occurs at the orbital occupied 
by the odd neutron. At fa ~ 0.3, the level crossing happens between 3/2[202] and l/2[330] 
orbitals, and leads to the discontinuity of adiabatic energy surface at fa m 0.3 in Fig. Ufa) 



as explained in the Ref. 30]. 



In Fig. [31 as an illustration, the neutron current in xy planes is plotted at z = 0.40 ~ 4.48 
fm for the ground state of 33 Mg. The direction and length of the arrows respectively represent 
the orientation and magnitude of current, which is in different scale for different panels. Take 
Fig. [3](a) as an example, since only azimuthal components of neutron current exist under 
the axial symmetry, the neutron current is on circular lines around the symmetry axis and 
peaks at r± = 3 ~ 4fm. From panel (a) to (f), one can see that, as z increases, the neutron 
current gradually grows larger and there is a tendency of closing up to the center in the xy 
plane. But in general, it peaks at r « R n = 3.4 fm (the neutron rms radius), indicating 
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FIG. 4: (Color online) The energy contributions due to the nuclear magnetic potential, Ej = 
—1/2 J d 3 r g U jU l p(r)j l Q(r), as a function of @2 for different configurations. 

that the current is mainly contributed from the old nucleon and flows around the nuclear 
surface 19] . It is interesting to note that for fixed z, the neutron current can be clock- wise, 
anti-clock-wise or even a mixture. 

In order to investigate the effect of the nuclear magnetic potential on the binding energy, 
in Fig. |4] the energy contributions for different configurations as a function of fli are given. 
In general, the nuclear magnetic potential makes the nucleus more tightly bound, and nu- 
merically, its contribution to energy is several hundred keV. Specially, the contributions for 

A, B, C, E are around -0.1 MeV, while for D, F, G -0.5 0.8 MeV, which is due to the 

breaking of proton pair shown in the following. 

The total energy (E cal ), quadrupole deformation parameter (/3 C ai)> valence nucleon con- 
figuration, the magnetic moments of the nucleus (^totai) and core (/i CO re) for different config- 
urations are listed in Table. [B It is obvious that the magnetic moments are sensitive to the 
valence nucleon configurations. For states A, C and E, the magnetic moments are negative, 
because the main component of corresponding valence nucleon wave function belongs to 
v I/7/2 which has a negative Schmidt value. For state B, the main component of valence 
nucleon wave function belongs to v and thus gives a positive magnetic moment. For 

states D and G, the total magnetic moments are remarkably large due to the breaking of 
proton pair and the enhanced Dirac current. For F, although there are three valence nucle- 
ons, the total magnetic moment is smaller than states D and G due to the strong cancelation 
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TABLE I: The energy (E ca \), the quadrupole deformation parameters (P C a,i), the valence nucleon 
configuration, the magnetic moments of the nucleus (^totai ) and the core (/i C orc)> in comparison 



with experimental value of the ground-state energy (E cxp ) 
The energy is in MeV and the magnetic moment in //n- 



321 ] and magnetic moment (/^exp) 



State 


E C al(E exp ) 


/3cal 


Valence nucleon configuration 


MtotaK/^exp) 


/^core 


A 


-251.85(-252.06) 


0.23 


A [330-] 


-0.9134(-0.7456) 


-0.042 


B 


-251.37 


0.40 


z^f [202+] 


1.6380 


-0.127 


C 


-250.69 


0.50 


4[321-] 


-1.6750 


-0.071 


D 


-250.58 


0.15 


z/i[330-] <g> 7r{|[211+]|[202+]} 


7.8179 


-0.216 


E 


-250.33 


-0.12 


z/f[303-] 


-2.0198 


-0.114 


F 


-250.18 


-0.11 


i>|[303-] <g> 7r{i[211+]|[211+]} 


1.4490 


-0.239 


G 


-249.34 


0.04 


i/|[312-] (8) 7r{|[211+]|[202+]} 


7.0082 


-0.288 



between the valence protons and neutron. 

The core polarization effect can be investigated by examining the magnetic moments of 
the core (/x cor e)- The baryon current of the valence nucleons is responsible for the polarization 
which results in the nonvanishing fi COTC . As it is seen from the Table. [U the fi core , which is 
caused by the polarization currents, is negative for all states. It is obvious that /^ C ore f° r 
states D, F and G are much larger than states ABC and E, due to the polarization effect of 
more valence nucleons. 



V. SUMMARY 

In summary, the configuration-fixed deformation constrained RMF approach with time- 
odd component has been applied to investigate the ground-state properties of 33 Mg with 
effective interaction PK1. Using the configuration- fixed deformation constrained calcula- 
tion, the ground state of 33 Mg has been found to be prolate deformed, fii = 0.23, with the 
odd neutron in 1/2 [330] orbital. The ground-state energy —251.85 MeV is close to the exper- 
imental value —252.06 MeV [32]. Using Eq. (Q, the magnetic moment — 0.913 /iN is obtained 
with the effective electromagnetic current which well reproduces the data — 0.7456 /in [ul 
self-consistently without introducing any parameter, in contrast with the shell-model re- 



10 



suits — 0.675 /xn and — 0.705 /xn [U| restricted to 2p — 2h configuration using two different 
interactions designed specifically for the island of inversion. The energy splittings of time 
reversal conjugate states, the nucleon current, the energy contribution, and the effect of core 
polarization due to the nuclear magnetic potential are discussed in detail. 

Apart from the core polarization, the meson exchange current correction is also very 
important for the descriptions of nuclear magnetic moment, especially the isovector magnetic 
moment. Investigations along this line are in progress. 
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